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FOREWORD 


This  report  deals  with  part  of  a  continuing  in-house  re¬ 
search  program  on  multi-component  flow  phenomena  conducted 
in  the  Energetics  Research  Laboratory  under  Project  No.  7116 
"Energy  Conversion  Research". 

The  author  wishes  to  express  appreciation  for  the  dis¬ 
cussion  he  has  had  with  Dr.  Hans  J.P.  von  Ohain,  Dr.  Brian 
Quinn,  and  Mr.  Siegfried  Hasinger.  Mr.  Brian  Hausfeld  is  to 
be  thanked  for  his  assistance  in  installing  and  operating  th 
experimental  apparatus. 
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ABSTRACT 


Velocity  of  sound  through  a  two-component  (solid-gas) 
medium  was  measured  as  a  function  of  solid  particle  concen¬ 
tration.  A  method  was  devised  in  which  phase  angle  between 
a  sound  source  signal  and  microphone  Bignal  could  be  measured 
as  a  function  of  solid  particle  density  in  the  apparatus. 
Acoustical  theory  provided  the  necessary  relationship  between 
phase  angle,  density,  and  sound  velocity.  Sound  velocity  in 
the  two-component  medium  was  measured  over  a  particle  volume 
fraction  range  from  10“^  to  0.5-  The  possibility  of  using 
this  technique  as  a  diagnostic  instrument  is  discussed. 
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INTRODUCTION 


Studies  of  sound  propagation  through  multi-component  media 
are  of  interest  for  a  variety  of  reasons,  A  better  understanding 
of  the  mechanics  of  multi-component  flow  can  be  obtained  from 
analysis  of  sound  propagation  through  the  multi-component 
medium,  both  from  direct  experimental  results  and  from  compari¬ 
son1  of  experimental  and  theoretical  predictions.  Experimental 
and  theoretical  results  indicate  that  sound  velocity  in  multi¬ 
component  media  is  dependent  on  particle  concentration,  which 
is  of  course  important  in  determining  Mach  r  imber  in  the  flow. 
Acoustic  instability  phenomena  are  of  direct  interest  to 

.1>.  -v  .  "  .  .  ..  .  .  .  ..  .  • 

feasibility  studies  of  the  colloid  core  reactor  concept 
(Ref.  .1).  Many  of  the  diagnostic  techniques  presently  employed 
in  the  study  of  multi-component  flow  are  at  a  rather  rudimen¬ 
tary,  level,  and  certain  acoustic  processes  offer  promise  as 
diagnostic  tools. 

A  considerable  amount  of  effort  has  been  devoted  to  the 
investigation  of  shock  waves  In  multi-component  media.  The 
first  analytical  work  in  this  urea  was  done  by  Carrier  (Ref. 

2),  and  Rudinger  extended  this  analysis  with  calculations  of 
shock  characteristics  in  two-phase  flows  (Ref.  3,  k ,  5). 
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The  results  showed  thnt  three  distinct  regions  exist  in 
the  shock:  t'ne  frozen  flow  zone,  where  gas-particle  transport 
processes  (drag  and  heat  transfer)  are  negligible;  the  relax¬ 
ation  zone,  where  gas  flow  and  particle  flow  interact;  and 
the  equilibrium  zone  where  normal  equilibrium  flow  is  again 
established.  In  the  frozen  flow  zone,  the  main  expansion 
shock  propagates  at  the  shock  velocity  in  clean  gas,  while 
the  particles  are  essentially  unaffected.  The  equilibrium 
wave  propagates  at.  a  velocity  that  is  dependent  on  the 
particle  concentration,  and.  consequently  the  Mach  number  of 
■the.  two-phase/,  flow  is  dependent  on  this  parameter.  Experi- 
mental  work  by  Eddington  ( R  f .  6)  has  confirmed  these  pre¬ 
dictions.  '  Using  water  droplet-air  mixtures;  hs  has,  shown 
that  shock  phenomena  -.occuring-  in  single  phase  supersonic  gas 
flow  is  duplicated  in  the  two-phase  flow,  with  appropriately 
modified  specific  heat  ratio  density,  dnd  sound  velocity. 

The  problem  of  sound  velocity  in  mi.ced-phds e  media-' 'has 
been  the  subject  of  a  number  of  investigations.  Sound  ve¬ 
locity  in  various  two-phase  media  have  been  calculated: 
phase  combinations  include  gas  bubbles  in  a  liquid  medium 
(Ref.  7),  solid  particles  in  a  liquid  medium  (Ref.  fl),  and 
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solid/liquid  particles  in  a  gaseous  medium  (Her.  5).  Jiira- 
scnic  resonance  has  been  need  t.n  verify  calculations  of  sonic 
speed  in  solid-liquid  mixtures  (Ref.  8,  9),  and  shock  studies 
have  verified  calculations  concerned  with  liquid-gas  flow 
(Ref.  6). 

The  experiment  described  in  this  report  is  a  determina¬ 
tion  of  the  sound  velocity  through  a  mixture  of  air  and  mod¬ 
erately  large  solid  particles  (glass  spheres  of  around  1+0  y 
diameter).  Maintaining  a  suspension  of  these  particles  in 
an  acoustic  resonance  tube  (the  technique  used  in  Ref.  8), 
for  sufficient  time  to  complete  the  measurement  is  difficult; 
therefore,  an  alternate  technique  is  used.  The  solid  particles 
are  allowed  to  fall  through  a  vertical  tube,  essentially  at 
terminal  velocity.  This  tube  is  part  of  an  acoustic  network 
containing  a  driving  source  (loud  speaker),  a  receiver 
(microphone)  and  an  acoustic  impedance  (connecting  tubes). 

Phase  difference  between  the  driving  signal  and  microphone 
signal  is  measured  as  a  function  of  particle  concentration 
and  the  relation  between  phase  difference,  particle  concen¬ 
tration  and  sound  velocity  through  the  powder  mixture  is  de¬ 
rived,  The  data  are  compared  with  calculated  values  of  sound 
velocity  as  a  function  of  particle  concentration. 
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concentration  suggest  a  possible  use  of  thiB  technique  In 
monitoring  particle  concentration  in  multi-component  flow. 
This  application  will  be  discussed  more  fully  in  Section 


II. 


DISCUSSION  OF  EXPERIMENT 


The  purpose  of  this  experiment  van  to  mens u re  the  speed 
of  sound  through  a  mixture  of  moderate  sized  particles  and 
air.  Solid  glass  spheres,  with  mean  diameter  around  UO  microns 
and  mass  density  of  2.35  gm/cm3,  were  chosen  for  the  experi¬ 
ment.  Particle  size  and  density  were  dictated  by  a  desire 
to  relate  the  results  of  the  measurement  to  investigations 
currently  being  conducted  by  this  laboratory  into  the  nature 
of  multi-component  vortex  flow.  The  resonance  technique  em¬ 
ployed  in  Ref.  8  is  probably  the  most  accurate  means  of 
measuring  sound  velocity.  That  method,  however,  was  not  well 
suited  for  this  experiment  because  of  the  difficulty  encountered 
-in  maintaining  a  stationary  suspension  of  such  large  particles. 

A  different  technique  was  devised,  and  is  explained  in  the 
following  text. 

Details  of  the  acoustical  system  are  shown  in  Fig.  1. 

Sound  vaves  generated  by  a  loud  speaker  are  transmitted  through 
the  tube  on  the  left  to  a  T  intersection.  A  portion  of  the 
sound  signal  continues  through  the  tube  on  the  right  and  is 
detected  by  a  microphone  placed  at  the  end  of  the  tube. 

The  vertical  tube  in  Fig.  1  acts  as  a  second  transmission 
line  connected  in  parallel  with  the  transmission  line  to  the 
microphone.  As  explained  in  Appendix  B,  an  acoustical  network 
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can  be  reduced  to  the  more  familiar  electrical  circuit  form 
by  definition  of  an  acoustical  impedance.  Fig.  2  and  3  show 
the  analogous  electrical  circuit  corresponding  to  the  acousti¬ 
cal  system  of  Fig.  1. 

Fig  2  is  the  complete  circuit,  showing  the  driving  voltage 
V  applied  to  the  loud  speaker,  the  transformation  from  electri¬ 
cal  energy  to  acoustic  energy  accomplished  by  the  loud  speaker 
(denoted  by  a  "black  box"?  )»  acoustic  impedances  of  the 
transmission  tubes  ( Z 0 ,  Zc,  and  Zp),  acoustic  impedance  of 
the  microphone  (  Zffl  ) ,  and  finally  the  conversion  from  acoustic 
to  electrical  energy  at  the  microphone  (t).  Under  certain 
conditions,  as  discussed  in  Appendix  B,  this  complete  circuit 
can  be  approximated  by  the  simpler  form  shown  in  Fig.  3.  The 
acoustic  impedance  of  the  vertical  tube  in  Fig.  1  is  denoted 
by  Zp,  acoustic  impedance  of  the  transmission  tube  to  the 
microphone  is  Zc,  and  the  acoustic  impedance  of  the  microphone 
is  assumed  to  be  purely  resistive  (F).  The  crux  of  the  experi¬ 
ment  is  the  measurement  of  phase  angle  between  the  A.  C. 
signal  driving  the  loud  speaker  and  the  signal  received  by 
the  microphone.  From  Fig.  3«  this  phase  angle  is  equivalent 
to  the  phase  angle  between  the  total  current  I  and  the 
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current  flowing  to  the  microphone,  Ig .  The  relationship  be¬ 
tween  phase  angle  and  the  circuit  elements  R,  Z  ,  and  Zp  is  de¬ 
rived  in  Appendix  B. 

Solid  particles  were  introduced  into  the  vertical  tube  .through 
the  funnel  shown  in  Fig.  1.  A  flow  diffuser  (a  cluster  of 
short,  small  diameter  tubes  placed  at  the  top  of  the  vertical 
tube)  provided  uniform  powder  flow  across  the  tube  cross-sec¬ 
tion.  The  glass  spheres  fell  through  the  tube  at  terminal 
velocity  (approximately  10  cm/sec),  and  uniform  density  was 
maintained  along  the  length  of  the  tube  with  constant  powder 
flow  rate  out  of  the  funnel.  As  shown  in  Appendix  B,  the 
impedance  Zp  is  dependent  on  the  density  of  the  air  column, 
p  ,  and  the  local  sound  velocity,  c.  By  measuring  particle 
concentration  simultaneously  with  the  phase  angle  between 
microphone  and  loud  speaker  signals,  the  sound  velocity  was 
determined  as  a  function  of  particle  concentration. 

A  dual  beam  oscilloscope  was  used  to  display  both  the 
microphone  signal  and  the  loud  speaker  signal.  The  speaker 
driving  signal,  supplied  by  a  pulse  generator  operating  at 
a  frequency  of  800  He,  triggered  both  beams  of  the  oscillo¬ 
scope,  so  that  the  phase  angle  between  driving  signal  and 
microphone  signal  could  be  read  directly  from  the  oscillo¬ 
scope.  The  phase  shift  produced  by  the  addition  of  solid 


7 


o  -?  r*  1  a  « 

r  —  *  -  - - ~ 


*  «  the 


CTlM  irv^i  jr  lU^aDui  XilQ 


the  displacement  of  the  microphone  signal  along  the  time  rxie. 
Large  magnifications  of  both  oscilloscope  axes  were  used  in 
order  to  facilitate  the  measurement  of  the  time  shift.  Pig. 

1*  shows  an  example  of  the  measurement:  the  peak  of  the  initial 
microphone  signal  (no  powder  flow)  is  superimposed  on  the 
shifted  signal  resulting  from  powder  flow  in  the  vertical 
tube.  The  time  increment  between  the  two  peaks  was  used  to 
determine  phase  shift. 

Particle  density  in  the  powder  stream  was  obtained  by 
two  different  methods:  light  absorption  for  low  density 
measurement,  and  mass  flow  rate  for  high  density  values. 

In  the  light  absorption  method,  the  powder  flow  was  positioned 
between  a  light  source  and  photometer;  the  decrease  in  light 
intensity  caused  by  the  flow  of  solid  particles  was  used  to 
calculate  the  number  density  of  particles.  For  high  flow 
density,  the  length  of  time  required  for  a  specified  mass 
of  powder  to  flow  through  the  tube  was  used,  along  with  the 
terminal  velocity,  to  determine  particle  concentration.  In 
the  overlap  region  (around  r.  *  .01),  both  methods  were 
used  to  insure  compatibility  of  the  two  sets  of  data.  Fig.  5 
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is  ft  plot  of  time  increment  LX  /ernus  particle  volume 
fraction  (  n  ) »  where  the  particle  volume  fraction  is  Just 
percentage  of  total  volume  occupied  hy  solid  particles,  and 
is  proportional  to  particle  mass  concentration. 
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III.  DATA  ANALYSIS 


*WC  U  a.  ti  t  1 t  i  6  S  V  £  i*  £•  Hi  c  d  £  u  i*  c  u  I  ii  t  li  1  o  cXpei'i- 


merit:  the  density  of  powder  flowing  in  the  powder  tube,  and 

the  phase  shift  of  the  microphone  signal  associated  with  the 
powder  density.  The  phase  shift  was  measured  on  an  oscillo¬ 
scope  by  observing  the  displacement  of  the  microphone  signal 
along  the  time  axis,  as  described  in  Section  II,  and  shown 
in  Fig.  1*.  Time  increment  is  related  to  change  in  phase 
angle  by 


At  =  (Mn  )  (1) 

2irv 

with  4)  denoting  the  new  phase  angle  between  speaker  signal 

and  microphone  signal,  is  the  phase  angle  with  no  added 

o 

powder,  and  y  is  frequency. 

The  relation  between  phase  angle  and  acoustic  parameters 
is  delineated  in  Appendix  B.  Specifically,  Eq.  (B21) 


tan  4> 


=  Lc^p/U3-LpR/oj 

Rp^/W^  +  Lp^+  R  RpTuP 


(2) 


I,c  ,  Lp  and  Rp  refer  to  the  circuit  elements  of  Fig.  3  and 
are  functions  of  acoustical  parameters:  R  must  be  determined 
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experimentally.  Lc  is  a  reasonably  simple  function,  Eq. 


(  B.  14)  : 


Lc 


cot 


U)i 

c 


(3) 


Since  p  ana  c  represent  conditions  in  the  free  air  of  the 
microphone  cavity  and  l  is  length  of  tube,  Lc  is  a  constant 
in  the  experiment.  With  the  appropriate  values  (£  =  25  cm, 
ui  =  5  x  1C)3  radian3/sec)  this  impedance  is 


L  =  -  .01  dyne-sec/cra3 

c  (M 

The  two  quant. ties  relating  to  powder  tube  impedance  are  more 
complicated.  As  derived  in  Appendix  B,  Eq.  (BlJ>) 


Rp  «  2ir2pc/S.y 


1  +  tan2  2tt£p/X 
1  +  UV4(a/j}‘H  tan*  2ir£p/X 


(5) 


a  1b  tube  radius,  A  is  wavelength,  and  £p 


With  the  values  a  =  2.5  cm,  A  =  1*0. 5cm  (corresponding 
•length  in  free  air),  £p  =  87  cm 


Rp/w  =  1.87  8p 


(6) 


with  8  ~  c /c ,  c0=  speed  of  sound  in  free  air. 
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The  Inductive  term  is 


_  On. 

03 


~  —  0  _  / 1 
ZL1  " 


I.  ~  0  /  A  / 


I  TTFTa7TTrt^ 


a/ _ u  an 


-?-'X 


21m  f 


(7) 


Upon  substitution  of  the  appropriate  values  in  Eq.  (7)j  again 
with  the  free  air  wave  length,  the  inductive  term  is 


bp  *  -8.95  Sp 

Eq.  (2)  can  now  be  written  as 


(8) 


tan  ^  * 


8.95  R/o)  -  .019 _ 

82.5  Bp  +  .089  ♦  1.87  R/o) 


(9) 


With  no  powder  in  the  powder  tube,  the  air  density  1b 
10_3gm/cm8,  and  3=1.  The  phase  angle  for  this  condition  was 
found  to  be  .9b  radian.  Putting  these  values  into  Eq.  (9) 
yields  the  value: 

R/o)  =  .  Ob  dyne-sec/cm^  (10) 

Finally,  the  phase  an^le  is  related  to  density  and  sound 
speed  by  the  function 


tan 


.  00b _ 

gp  +  .002 


(11) 
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Figure.  5  shows  experimental  measurements  of  time  shift 
At  as  a  function  of  powder  density.  Dots  on  this  curve  re¬ 


present  density  values  obtained  from  the  light  absorption 
technique,  while  crosses  refer  to  data  from  the  mass  flow 
rate  method.  (For  a  complete  discussion  of  these  techniques, 
refer  to  Section  II).  The  solid  curve  is  calculated  from 
Eq.  (11),  assuming  8=1.  The  calculated  curve  approaches 
experimental  results  in  the  limits  of  low  volume  fraction 
and  high  volume  fraction,  but  between  those  limits  a  large 


discrepancy  exists.  Assuming  this  discrepancy  is  due  to 


the  function  8(  h),  Eq.  ( 1 1 )  can  be  applied  to  the  experi¬ 
mental  points  in  Fig.  5  to  calculate  sound  velocity.  This 
result  is  Bhown  in  Fig.  6,  with  the  theoretical  function  de¬ 
rived  in  Appendix  A  shown  as  a  solid  o.urve. 

Agreement  between  experimental  and  theoretical  sound  ve¬ 
locity  1b  fairly  good.  The  random  variations  that  are 


apparent  in  Fig.  6  are  due  mainly  to  experimental  errors  in 
measuring  powder  density,  as  errors  involved  in  measuring 


phase  shift  should  be  negligible.  Some  of  the  errors  inherent 
in  the  density  measurements  are  as  follows:  non-uniform 
density  through  the  powder  tube,  irregularities  in  powder 
feed  (which  would  give  variation  of  instantaneous  powder 
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density  that  would  not  be  seen  from  integrated  powder  flow 
rate  measurement)  ,  and  the  slow  response  time  of  the  photo¬ 
multiplier  (of  the  order  of  .5  Bee)  which  would  mask  rapid 
fluctuations  in  powder  density.  There  is  a  better  agreement 
between  theoretical  values  and  data  taken  with  flow  rate 
determination  of  density  because  each  of  these  points  is 
actually  an  average  of  5  independent  measurements:  thus 
most  of  the  random  spread  has  been  removed  from  these  data. 

The  point  at  =  .4  shows  a  significantly  large  deviation 
from  the  calculated  curve.  At  such  large  particle  concentra¬ 
tion,  particle-particle  interaction  may  be  of  importance, 
noticeably  altering  the  terminal  settling  velocity  for  the 
particles.  Also,  calculation  of  ^becomes  quite  sensitive 
to  error  in  measuring  time  shift.  At,  as  well  as  initial 
phase  angle,  at  high  particle  concentration.  Therefore, 
systematic  error  may  be  introduced  in  this  region,  and  less 
confidence  should  be  placed  on  the  data  at  high  concentrations. 

A  systematic  trend  below  the  theoretical  curve  seems  to 
be  in  evidence  in  the  vicinity  of  a  volume  fraction  of  .005* 
This  trend  may  be  the  result  of  an  approximation  made  in 
arriving  at  Eq.  ( 11  )  .  The  constants  in  Eq.  (11)  are  really 
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combinations  of  Lp  (from  Eq.  (7)  and  (8))  and  Rp  (Eq.  (5) 
and  (6)).  The  approximation  involved  is  the  neglect  of  the 
dependence  of  Lp  and  Rp  on  wavelength.  The  complex  depen¬ 
dence  makes  a  complete  analysis  of  these  parameters  both 
complicated  and  sensitive  to  errors  in  8.  The  sharp 
variations  in  Rp  and  Lp  with  wavelength  occur  due  to  resonances 
between  an  incoming  sound  wave  and  a  wave  reflected  from  the 
open  end  of  the  tube.  The  introduction  of  powder  into  the 
tube  may  smooth  out  these  sharp  resonances,  as  each  particle 
can  act  as  a  source  of  reflected  sound.  Then  some  averaged 
value  of  Rp  and  Lp  would  be  appropriate.  This  effect  may 
explain  the  surprisingly  good  agreement  between  theory  and 
experiment  obtained  by  completely  neglecting  wavelength 
variation  in  the  calculation  of  the  parameters  R  and  L  . 

r  Jr 
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IV.  CONCLUSIONS 


Results  of  this  experiment  indicate  that  the  propagation 
velocity  of  sound  in  a  solid-gas  medium  is  dependent  on  the 
concentration  of  solids,  and  the  sound  velocity  for  highly 
loaded  mixtures  may  he  only  a  small  fraction  of  the  pure  gas 
sonic  velocity.  The  theoretical  calculation  of  velocity 
ratio,  8,  in  Fig.  6,  shows  a  minimum  value  of  . OU  at  a 
particle  volume  fraction  of  0.5*  while  experimental  data 
reaches  a  minimum  value  at  a  slightly  lower  volume  fraction. 
Beyond  this  point,  8  increases  rapidly  with  increasing 
particle  load.  According  to  Eq.  (A35),8  becomes  infinite 
at  e«o,  but  in  reality  this  limit  should  be  the  sound  velocity 
in  the  solid  material,  around  50JO  m/sec  for  glass.  This 
minor  discrepancy  is  due  to  the  neglect  of  the  compressibility 
of  the  solid  particle.  The  higher  relative  velocity  obtained 
experimentally  at  U=  may  be  the  result  of  sound  trans¬ 
mission  through  clusters  of  solid  particles.  At . such  high 
particle  concentration,  local  irregulariti'.-c'  ’.’h  fixture 
density  may  approach  the  solid  material  den.  tv  .  and  sound 
transmission  through  this  region  would  proceed  at  a  much 
higher  rate  than  through  the  bulk  distribution.  No  con¬ 
clusive  statement  can  be  made  in  this  regard  at  present. 
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v.  o  r*  Vifr-onaia  nf  t.hf»  ti  n  r.  p  r  t.  a  1  n  t  v  Involved  in  this  datim 
point  (refer  to  the  comments  in  Section  III). 

In  deriving  the  expression  for  sound  velocity  in  the 
solid-gas  mixture,  it  was  tacitly  assumed  that  particle 
inertial  effects  were  negligible;  i.e.,  it  was  assumed  that 
solid  particles  were  moved  by  the  pressure  wave  identically 
with  the  gas,  so  that  density  of  the  medium  could  be  taken 
as  the  sum  of  densities  of  the  two  components.  The  general 
agreement  between  experimental  data  and  the  theoretical 
curve  seems  to  justify  this  assumption,  but  with  particles 
of  larger  diameter,  or  higher  density,  inertial  effects 
may  become  important. 

Analysis  of  many  multi-component  flow  phenomena  would 
be  significantly  simplified  if  one  could  consider  the  in¬ 
homogeneous  medium  as  a  perfect  gas,  with  some  modification 
of  the  appropriate  parameters.  As  an  example  of  this  simpli¬ 
fication,  Eq .  (A20)  can  be  used  to  obtain  an  expression  for  3 
in  the  solid-gas  mixture.  The  modified  specific  heat  ratio  r 
is  substituted  for  Y«  and  the  density  of  the  mixture  is  taken 

asco  +  (l-e)n  •  The  result  for  the  ratio  of  sound  velocity 

g  P 

in  the  mixture  to  that  in  the  perfect  gas  is 
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LorapariBon  oi  uus  resuit  vun  me  more  exacting  aerivation 

in  Eq.  (A35)  shows  that  the  two  results  differ  only  by  a 
-1/2 

factor  e  .  These  two  functions  are  plotted  in  Pig.  8: 

Eq  (12)  is  shown  &b  a  solid  curve,  and  Eq.(A35)  is  shown  with 
a  broken  curve.  Since  e  *lfor  small  values  ofh(n«  1-E)*  the 
two  curves  are  practically  equivalent  up  to  a  solid  volume 
fraction  of  0.01.  Thus  it  appears  that  these  simplifying 
assumptions  can  be  used  safely  for  particle  volume  fraction 
lower  than  about  .01. 

Vortex  flows  with  sub-micron  size  condensation  particles 
have  been  studied  by  Fletcher ,  et  al  (Pef.  13).  Volume 
fraction  of  liquid  in  these  experiments  was  on  the  order  of 
10~\  so  that  sound  velocity  should  be  around  90 of  the 
value  for  pure  gas,  and  Mach  number  of  the  flow  would  be 
about  10^  higher  than  the  corresponding  value  for  pure  gas 
flow.  Investigations  into  the  nature  of  vortex  flows  with 
high  solid  particle  loading  are  currently  being  conducted 
in  this  laboratory.  Solid  volume  fraction  aB  high  as  0.1 
has  bren  obtained  in  this  work,  with  tangential  veloc it ies 
considerably  higher  than  the  20  m/sec  sonic  velocity  pre¬ 
dicted  by  Fig.  6.  Study  of  these  flow  conditions  with  high 
loading  may  provide  further  insight  into  the  dependence  of 
sound  velocity  on  particle  concentration. 
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With  only  minor  modification!  the  experimental  apparatus 
described  in  this  report  could  be  used  as  a  diagnostic  in¬ 
strument  for  measuring  density  in  linear  flows.  The  powder 
flow  to  be  measured  would  be  directed  through  the  vertical 
tube  of  Fig.  1,  and  the  resulting  phase  shift  converted  to 
particle  density  (or  volume  fraction)  with  the  calibration 
curve  in  Fig.  5«  One  problem  associated  with  this  technique 
may  be  interference  from  noise  generated  by  the  powder  flow. 
More  Btudy  must  be  given  this  point  to  find  a  suitable  opera¬ 
ting  frequency  with  little  extraneous  noise,  or  develop  a 
noise  filtering  system.  Of  particular  interest  to  this 
laboratory  is  the  development  of  diagnostic  equipment  suit¬ 
able  for  use  with  swirl  chambers  such  as  that  described 
in  Ref.  13.  Basically,  these  chambers  are  large  diameter 
cylinders  with  axial  vortex  flows.  Boise  level  during  opera¬ 
tion  of  these  chambers  is  rather  high,  so  noise  interference 
would  be  a  major  problem  for  the  phase  shift  measurement. 

The  particle  concentration  obtained  from  the  phase  shift 
measurement  would  probably  be  an  average  over  the  entire 
vortex  flow,  as  the  whole  chamber  would  be  encorporated  in 
the  acoustic  system.  However,  Ref.  lU  indicates  that  at 
sufficiently  high  flow  rate,  sound  attenuation  perpendicular 
to  the  flow  becomes  strong.  Under  these  conditions,  localized 
density  measurement  with  the  phase  shift  experiment  may  be 
possible . 
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VI. 


APPENDIX  A 


„  t  C„ _ A  C.. - J- J 

~t ~  ...  W^B..  ~  ~  ^ 

In  this  section  the  equation  of  propagation  of  a  pressure  wave  through 
a  suspended  powder  is  developed.  In  order  to  compare  the  result  directly 
with  the  pure  gas  solution,  the  classical  derivation  for  a  homogeneous  gas 
is  first  presented  (taken  largely  from  the  derivation  by  Morse,  -<ef.  10). 

We  consider  the  propagation  of  a  plane  pressure  wave  through  a 
cylindrical  tube  of  cross-sectional  area  S.  |(x,t)  denotes  displacement  of 
gas  due  to  the  wave,  P  is  pressure,  and  P  is  gas  density. 

The  total  mass  within  the  volume  Sdx  is 


dm  ■  pSdx 


<A1.) 


This  quantity  is  conserved  as  the  pressure  wave  passes.  As  the  waves 
pass,  one  bounding  plane  moves  a  distance  £  (x),  while  the  other  plane 
moves  a  dista#C*C(x+  dx).  The  volume  element  changes  to  a  value 


Sfdx-^x)  ♦  £(x+dxl] 


Sdx  +  s|j^x 


(A2) 


The  density  is  taken  to  be 

p(x,t)  ■  p0 ( 1*6 ) 


(A3) 


where  6  is  considered  to  be  much  less  than  1.  With  conservation  of  the 
mass  in  Eq,  (Al)  we  have 


dm  ■  P0Sdx  ■  SP0(l+fi) 


(A4) 
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This  equation  gives  the  resul..  o  first  order,  that 


The  equation  of  state  for  a  perfect  gas  is 
PV  ■  RT 

A  temperature  differential  d  T  gives 
RdT  ■  PdV  +  VdP 

and  dividing  Eq.  (A7)  by  (A 6)  yields 

dT  dP  dV 
T  *  V  *  V 


(A5) 


(A6) 


(A7) 


(A8) 


Heat  capacity  at  constant  volume, C^,  is  defined  as  the  amount  of  heat, 
dQ,  absorbed  per  unit  temperature  increment  with  constant  volume 


.  I ®/dP\ 

\dT/v  constant  3rldTiv  constant  (A9) 


,  so  we  have 

30  F 
3P  T 


(A10) 


In  a  similar  manner  we  derive  the  heat  capacity  at  constant  pressure: 


c  .  ia  v 

P  3  V  T 


(All) 
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2k 


With  mass  PSdx  and  acceleration 


37? 


,  the  equation  of  motion  is 


"H**  '  If 


(Aia) 

Combining  Equations  (A18),  (A16)  and  (A5)  we  have  the  two  wave  equations 


YPo 

at*  -ifr  57? 

-So 

p  ax* 

These  equations  indicate  ^hat  pressure  and  particle  velocity  are 


propagated  by  wave  motion  with  the  phase  velocity 

CL 


l -Ho 
°o 


( A19) 


( A20  ) 


The  result  in  Eq.  (A20)  is  the  sonic  speed  in  a  homogeneous  gas: 

now  we  shall  revise  the  derivation  by  assuming  that  solid  particles  are  suspended 

in  the  gas,  with  void  fractions  being  defined  as  the  ratio  of  gas  volume  to 

total  volume,  and  particle  volume  fraction  n  defined  as  the  ratio  of  particle 

volume  to  total  volume.  The  ratio  of  particle  density  P  to  gas  density  p 
•  P  * 

is  *m  * 

The  total  mass  in  the  test  volume  Sdx  is  now  given  by 
dm  «  Sdx[ (l-e)pg  +  epp] 

am  .  r„li-e  Jjji]pss4x  (A21) 

The  equation  of  state  for  the  solid-gas  mixture  will  be  written  as 


P(V-Vp) 


RT 


(  A2d) 


with  the  result  corresponding  to  Eq.  (B8)  of 

.  dV  .  dT 
f  7-Vp  “? 

or  dP  dV  dT 

7  +  ev  ■  “? 


(  A23) 


The  heat  capacity  calculations  are  changed  somewhat  with  the  addition 
of  powder,  since  the  solid  heat  capacity  Cg  must  now  be  considered. 
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t\l  constant  voiumr,  me  amuum  ui  neat  auueu  pet  umc  rise;  u  icuipctatutc 


(dfi)v-(Cw*C„)dT 


(A2U) 


and  from  Eq.  (423), 


Similarly, 


ft  t3v*=v*C, 
ft  “  (cv  +  c8)  _1 

(dQ)p-(C  ♦C>)dT 


(A25) 


which  leads  to  the  result 


■  “v=.> 


(A26) 


When  both  P  and  V  vary,  the  heat  addition  is 

48  ■ .  ft  4,r  *  ft" 

40  -  [  ( Cp  ♦  c,)  |1  *  (cv  ♦  o.)  <AST> 

and  the  adiabatic  restriction  yields 

H  U86) 

and  with  the  result  of  Eq;  (A2),  we  have 

-  -cOt+C,  P 

3*  Cp*C8  P0 

_  -r  lJt+Cm/Cm  P 

lx  l+Cg/Cp  Pe  (A29) 

If  hp  s  heat  capacity  of  gas  at  constant  pTessure,  and  h8  :  heat  capacity 


(A28) 


of  solid  CB  ■  Pp  (1-C)hs 

Cp  ■  0g  Ghp 

B  =  Cs/Cp  ■  r„  i^£  jli. 

nP 

and  Eq.  (A29)  can  be  expressed  in  the  simpler  form 

*.,1/112  _L 

ax  *  l+B  PQ 

P  -  £>od 

with  1»B 

r  17*51” 


(A30) 


( A31 ) 
(A32) 
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The  equatiop  of  motion,  analogous  to  Eq.  (A18)  is 

h  &  *  -  -If 

L 


(A33) 


and  we  have  the  wave  equation 

325«  rp» _ 

"Tt*  ermPg|il-erm-l 

r*  . 


(A3M 


1  rpo 

i  ""1 

J  Crmp 

L.tJtel  Ogl 
f  ,nm  J 

c/c0- 

Jm  > 

V  ^ejl-  e  £,-1 

V  ferTTfr.] 


f  /  Y*  1  »  B 
1  +  BY 

For  the  particular  c  ase  of  solid  glass  spheres  suspended  in  air,  we 
have  the  following  parameters: 
rra  -  1.9  x  103 


y  -  i.k 

h,  •  0.12 
hp  ■  0.22 


(A35) 

(A36) 

(A3T) 


( A  3  Q  ) 


B  ■  1  x  10 


i  l-e 


( A39 ) 


B/y  -  715 


The  function  (3  is  plotted  in  Fig.  6,  and  T  is  plotted  in  Fig.  7. 
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VII.  APPENDIX  B 


Acoustic  Impedance  and  the  Equi valent  Circuit 


In  electrical  circuit  theory,  a  quantity  called  the  electrical  impedance 
of  a  circuit  element  is  defined  as  the  ratio  of  applied  voltage  to  current, 
and  it  may  be  a  complex  number.  Well  known  rules  exist  for  calculating  the 
effect  of  combinations  of  individual  circuit  elements  (e.  g.  ,  series -parallel 
combinations  of  resistors,  capacitors,  or  coils)  so  that  electrical 
characteristics  of  complicated  circuits  can  be  deduced  by  proper  application 
of  those  rules.  In  this  appendix,  the  acoustical  impedance  is  derived, 
enabling  electrical  circuit  theory  to  be  applied  by  analogy  to  acoustical 
network  problems.  The  derivation  is  essentially  taken  from  Morse,  (Ref.  10) 
The  pressure  p  in  the  acoustical  system  will  correspond  to  electrical 
voltage,  particle  velocity  5  corresponds  to  electrical  current,  and 

the  acoustical  impedance  is  therefore 

Z  -  \ 


Consider  propagation  of  a  sound  wave  through  a  tube  of  radius  a  and  length  1. 
As  shown  in  Appendix  A,  pressure  and  particle  displacement  can  both  be 
expressed  as  wave  functions.  For  a  tube  of  finite  length,  a  reflected 
wave  as  well  as  the  primary  wave  will  be  present,  so  that  pressure  function 


will  have  the  form 


p  -  P+.  ik(x-ct>  ♦  Vm9  -**(x+et) 


P+  is  the  magnitude  of  the  primary  wave,  ]P.  is  magnitude  of  reflected 
wave,  c  is  the  propagation  velocity,  k  *  w/d,  and  w  is  propagation 
angular  frequency.  The  complex  number  ♦is  defined  such  that 

P-/P+  •s'2"  (V 


2b 


Then  p(x)  ■  P 


-♦-iwt  (  kx  >M>lkx\ 

*•  v*  .  I 

p(x)  ■  2P+#’lf“1W1i  sinh 


In  a  similar  fashion,  particle  displacement  can  be  written  as 


i(x )  ■  P 


fc+ikx  A.-W-ikx) 


where  the  constant  K  must  be  determined. 

From  Eq.  (Al6)  and  (A20)  (or  equivalently  Eq.  (A31)  and  (A35)  pressure 
and  displacement  are  related  by 

H 

Evaluation  of  this  equation  at  x»0  gives  the  complex  constant  value 

K  •  (b6) 

pew  v 

Particle  velocity  is 

t  m  i  (alt»+ikx  +a-V-ikx)  (»tw) 

Pew 

t  ■  i£k  .  -^imi  <*♦*»*> 

PC  c  \ by  I 

Therefore,  in  accordance  with  Eq.  (Bl) 

Z  ■  Po  tanh  (U  +  (B8) 

This  equation  will  now  be  applied  to  two  specific  problems  of  interest, 
namely  a  calculation  of  acoustic  impedance  of  a  closed  tube  of  length  1, 
and  an  open  tuba  of  length  1, 

As  mentioned  previously,  ♦may  be  complex.  Therefore,  let  us 

define  the  quantity  y  «.  _&££,  ■  w(a-ib  }  (B9) 

c 

a  •  -lReb>  ,  b  -  — jj.In  ♦  ■ 


The  ratio  of  the  magnitude  of  reflected  wave  to  incident  wave  is 
Now  suppose  that  the  propagation  tube  is  closed  at  x*l  by  a  rigid  plate 
Then  the  magnitudes  of  reflected  andincident  wave  are  equal,  giving  a»a 
The  impedance  term  reduces  to 
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7  ■  Pc  t*nh(-i1Tb)  *  -iPc  tanitb 


(BIO) 


a.  1  *  U  A  J  _  J-l _ 

“  A  I  y  Mb  *tuy  v  Uttlt'a  W  A  O  AMI  I  A  «^Uli  AM|( 


1/2  ■  -  JL 


wl 

tic 


or 


-I_V 

m 


1  ^  Si 

2  c 


(Bll) 


Then 


1/2  ♦  Si.  _K* 

'  irc  ""c 


The  impedance  of  the  closed  tube  at  x  «o  is  therefore 


Z0  »  -iPc  ta 


»  (*♦■*) 


or 


z  ■  i Dc  cot  Si. 

O  A 


( B12  ) 


The  derivation  of  impedance  of  an  open  tube  is  too  lengthy  for 
presentation  here:  reference  should  be  made  to  Ref.  10  for  details. 

The  acoustic  impedance  at  the  driving  snd,  x*o ,  of  an  open  tube  of 
length  1  and  radius  a  is  given  approximately  by 

z  «-ipc  g.ng./**  igy.ls/i)*  „ 

°  1  -  tan  2*1  /X  (B13) 


with 


lp  ■  I  ♦  8a/3ll 


The  acoustic  impedances  derived  above  can  now  be  applied  to  the  acoustical 
system  shown  in  Fig.  1.  As  outlined  in  Ref.  11,  an  equivalent  electrical 
circuit  can  be  constructed,  and  is  shown  in  Fig.  2.  Coupling  between 
electrical  and  acoustical  energy,  as  performed  by  both  the  loud  speaker 
and  microphone,  is  denoted  by  T;  Zmis  mechanical  impedance  of  microphone, 
and  will  be  assumed  purely  resistive.  ZG  represents  impedance  of  the  tube 
connecting  the  loud  speaker  to  the  T  section,  Ze  is  impedance  of  tube  from 
T  section  to  microphone,  and  Zp  is  impedance  of  powder  tube.  The  frequency 
was  selected  so  that  the  phase  shift  resulting  from  Z0  was  II  radians,  and 
assuming  equivalent  reactance  in  coupling  at  the  loud  speaker  and  microphone, 
the  equivalent  circuit  in  Fig.  2  can  be  reduced  to  the  circuit  in  Fig.  3, 
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with  the  stipulation  that  a  phase  lag  of  JI  radians  is  introduced  in  the 
signal  generaiux . 

From  Eq.  B12,  the  impedance  Zc  (closed  tube)  is  purely  inductive,  i.  e. 


i*»lc 

£C  cot«i 


(B1U) 


Eq.  B13  indicates  that  Bp  (o^en  tubi)  is  complicated:  after  some 

algebraic  manipulation  the  impedance  can  be  expressed  in  the  simpler 

form  ZB  ■  R 

p  P  P  a 

Rp  -  2nVUA)J  - rs-T-n  (B15) 

p  1  +  4 tt  (a/J)  tan  2tt l_ / X 


-  _££ 
¥ 


tan  gitlp/X  »  Uirl>(a/X)l,tan  2trlpyx 
1  ♦  Uir1*  ( a/*  J^tan^iTip/). 


Currents  in  the  parallel  circuits  are  divided  such  that 

Ij_  Zp  •«  i2(a  +  ze)  { Bl6 ) 

'll  (H  +  lwLp)  ■  I2(R  +  ivLe ) 

(I  -  Ig)(Pp  +  lwLp )  »  I2(R  ♦  ivLc)  (BIT) 


Rearrangement  of  terms  leads  to 

■=£■  ■  Rp  ♦  IvLp 

R  +  Rp  +  i»»(Lc  +  Lp) 


(318) 


By  multiplying  both  numerator  and  denominator  by  the  complex  conjugate  of 
the  denominator,  the  complex  term  in  the  denominator  is  removed,  with 


the  result 

~  -  L( R*Rp)Rp»»a Lp ( LC4LP ) ]  *  i»[ Lb ( R*Fp ) -Rp ( Lr  +  l* } 
1  (R+Rp)*  w2  (  Lc  +  Lp  )  2 
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^  and  Its  obtained  from  the 

(B20) 


where  of  course  Im  denotes  the  imaginary  portion  of  the  function  and 
Re  the  real  part. 

tan  $  ■  -  »[Lp(R  +  Rp)  -  Rp(Lr  +  LPl] _ 

( R  +  Rp )Rp  +  w2  Lp( Lc  +  Lp ) 

or 

tan*  -  ka.Rp/w  r.-kpA/JS -  (B2i) 

Rp2/w2  ♦  Lp2  +  P  Rp/u2  +  LpLc 

The  phase  angle  represents  the  phase  lag  between  currentsl|  and  I,  or 
equivalently  the  phase  lag  between  speaker  signal  and  microphone  signal. 
Therefore  Eq.  (B21)  is  the  relationship  between  the  acoustical  properties 
of  the  system  and  the  measured  phase  angle  required  to  complete  the 
analysis  of  the  experiment. 

The  functional  forms  of  Lp,  Rp,  and  Lc  have  been  derived:  only  R 
remains  to  be  determined.  Since  R  is  the  microphone  resistance  and  will 
contain  acoustic  as  well  as  electrical  terms,  it  is  best  determined 
experimentally,  as  is  explained  in  the  data  analysis  section. 

'This  formula  differs  from  the  standard  relation  given  in  texts 
on  electrical  circuit  theory  by  the  factor  -1  (see  for  instance 
Ref.  12).  The  minus  sign  is  due  to  a  difference  in  writing  the 
oscillating  time  dependance,  which  for  standard  electrical  circuit 
theory  is  written  as  ,  as  compared  to  the  form  e-iwt  used 

in  this  derivation.  The  result  is  a  difference  in  sign  for  the  phase 
lag  $  . 


The  phase  angle  between  currents 


equation 


tan<$  ■  -  Im (  1 2 / 1  ) 
Re(  I.,/t  1 
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hift  as  Function  of  ^article  Volune  Tractior 


Figure  7.  Ratio  of  Specific  Heats  as  Function  of  ’’’article  Volure  Fraction 


—  EXACT  CALCULATION 

—  NEGLECTING  PARTICLE  VOLUME 


of  ^article  Volume  Fraction 
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